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Abstract 

Motivated by the unique physical properties of biological active matter, e.g., cytoskeletal dy- 
namics in eukaryotic cells, we set up effective two-dimensional {2d) coarse-grained hydrodynamic 
equations for the dynamics of thin active gels with polar or nematic symmetries. We use the well- 
known three-dimensional {3d) descriptions [K. Kruse et al, Eur. Phys. J E, 16, 5 (2005); A. Basu 
et al, Eur. Phys. J E, 27, 149 (2008)] for thin active gel samples confined between parallel plates 
with appropriate boundary conditions to derive the effective 2d constitutive relations between ap- 
propriate thermodynamic fluxes and generalised forces for small deviations from equilibrium. We 
consider three distinct cases, characterised by spatial symmetries and boundary conditions, and 
show how such considerations dictate the structure of the constitutive relations. We use these to 
study the linear instabilities, calculate the correlation functions and the diffusion constant of a 
small tagged particle, and elucidate their dependences on the activity or nonequilibrium drive. 
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I. INTRODUCTION 



The dynamics of equilibrium systems describe time evolutions of fluctuations around 
the minimum free energy or maximum entropy states. Systems in equilibrium do not con- 
sume energy continuously. Dynamics in such systems are characterised by the Fluctuation- 
Dissipation-Theorem (FDT), relating a susceptibility with an appropriate correlation func- 
tion In contrast, active systems are driven out of equilibrium by a continuous consunip- 
tion (supply) of energy. Well-known examples are gels driven by chemical reactions |2|, 
vibrating granular materials jsl, large scale coordinated motion of self-propelled particles 
like bacteria colonies and bird flocks ^ and cell cytoskeleton [5]. The physically interesting 
aspect of biologically relevant active matter is its ability to convert free energy (available in 
the form of chemical energy, e.g., hydrolysing Adenosine- Triphosphate (ATP)) into mechan- 
ical work and systematic movement. Many of these systems, despite having very different 
characteristic length and time scales, and microscopic details, share common general features 
in the long wavelength, large time limit where conservation laws and symmetries (and not mi- 
croscopic details) govern the general behaviour. Such generalities and hence coarse-grained 
approaches based on them are therefore generic. Coarse-grained descriptions in statistical 
)hysics have a long and successful history, beginning with equilibrium critical phenomena 
6| and equilibrium critical dynamics and more recently, to a variety of systems out of 
equihbrium js, 9|. Of late, such coarse-grained approach has been applied extensively to 
understand the physics of biologically motivated systems, see, e.g., instabilities of cortical 
actin layer 10|]. Nonequilibrium nature of the fluctuations in cell, and hence violation of 
the FDT, has been tested experimentally by Mizuno et al [11]. Hydrodynamic theories have 
been developed which successfully describes some dynamic cellular processes 12|. Further, 
formation of specific patterns in cytoskeletal structures is common, e.g., bundles, asters, 
vortices etc. [13|. Ref. IJ] provides a coarse-grained theory at the mesoscopic scale for 
cortical patterns in plant cells. For further details, we refer the reader Refs. 15l-ll7l| for 
recent reviews on these subjects. 

In this paper we focus on the dynamical behaviour of a cortical layer of filaments un- 
der various circumstances (e.g., various spatial symmetries, boundary conditions) at length 
scales much larger than the filament lengths and layer thickness with polar as well as ne- 
matic macroscopic ordering, for which a generic coarse-grained continuum 2d description 



would be appropriate. Our main achievement is a set of coarse-grained hydrodynamic equa- 
tions of motion for the orientational order parameter, local density and the velocity fields 
in the frictional limit, which we use to obtain results on the macroscopic properties of the 
underlying systems subject to various conditions. We consider the three following distinct 
cases for a thin active gel: (i) In-plane nematic or polar order without any external force 
(hereafter System I, (ii) No in-plane order (ordering is normal to the plane) (hereafter 
System II and (iii) In-plane polar order with external (surface) forces (hereafter System 
III). In each of the cases, we set up the effective 2d equations by using the 3d framework of 



Refs. 



18 



19| for 3d active gels and averaging them over the thin direction. We show that 



such effective 2d equations are consistent with the macroscopic symmetries of the systems. 
We show that our equations for finite activity (characertised by a constant parameter A/x 
in our notation) break time-reversal invariance explicitly and hence display nonequilibrium 
behaviour. In each of the the magnitude of Afi < exceeds a critical value A/i^, 

the initial chosen states exhibit finite wavevector instability, leading to anisotropic patterned 
states. Below the threshold of instability, i.e., |A/i| < |AyU,c|, equal-time density correlations 
function display giant fluctuations in some of the cases considered by us, as also reported 
previously in Refs. js, |2o| for active nematics and polar flocks. In addition, we flnd generic 
underdamped propagating modes for systems with in-plane polar order both below and 
above the threshold of linear instability. Our fluctuating hydrodynamic theory of 2d active 
matters in contact with a substrate should be well-suited to describe a thin actin cytoskele- 



ton in contact with the bulk of the cell, which serves as a substrate. Refs. [18|, |21[ briefly 
discussed dynamical properties of a quasi-2(i active gel layer in contact with a substrate. 
Here, we extend their work and provide systematic studies of quasi-2(i thin active gel layers 
in contact with substrates together with various boundary conditions and initial reference 



states. We show how effective friction, which dominates the ve^ 



ocity fleld dynamics, emerges 



in each of the cases out of the lubrication approximation [22| we make here. In addition, 
the two-dimensional description of actin dynamics in a motility assay is related to one of 
the cases considered by us here. The problems which we address here broadly refer to the 



phenomenon of flocking: the collective motion of self-propelled active particles; see Ref. 



23| 



for a review of recent results. We use the framework developed in Ref. [18| (equivalently, 
Ref. jsl). The remaining part of this paper is structured as follows: In Sec. |TT]we discuss 
the appropriate 2d thermodynamic fluxes and conjugate generalised forces for this problem. 
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In Sec. Illlt we derived the constitutive relations between them. Then in Sec. [IV] we add 
FDT-obeying thermal noises in the 2d equations when A/x = 0. In Sec. |Vl we use our consti- 
tutive relations to illustrate linear instabilities at finite wavevector and calculate correlation 
functions of the local polarisation and density variables. We further calculate the diffusion 
coefficient of a tagged particle and highlight its non-trivial dependence on the activity. In 
Sec. I VI I we summarise and conclude. Finally, a short comment on the notations used in the 
article is in order: We use a Roman subscript (e.g., i) to denote components of a 2d vector, 
i.e., i = x,y and a Greek subscript (e.g., a) to denote the components of a 3d vector, i.e., 
a = x,y, z. 



II. GENARALISED FLUXES AND FORCES 

We develop a dynamical theory of active gels in the coarse-grained (continuum) hydro- 
dynamic approach. We consider here a single fluid model where the embedding solvent is at 
rest [19|. The system is assumed to be only slightly away from equilibrium. This available 
free energy due to the active, nonequilibrium processes is equivalent to a nonequilibrium 
generalised force, which breaks the time reversal symmetry, represented by A/i in our de- 
scription. The linear constitutive relations between the appropriate fluxes and generalised 
forces, which we obtain as a general expansion of the fluxes in terms of the forces in the spirit 
of Onsager reciprocity relations for fluctuations in equilibrium, holds for active gels close to 
thermodynamic equilibrium. The detail form of the linear expansion of the fluxes in terms 
of the forces depend upon the symmetries of the system under consideration (see below). 
We consider the dynamics in the high friction limit. Hence, the constitutive equations we 
set up are valid only in a preferred frame of reference and as a result there is no Galilean 
invariance of the system - they will not be invariant under the transformation v — )■ v + vq. 
Here v is the local velocity of actin gel filaments. Further since we are modeling a gel with 
local nematic or polar order, we introduce a polarisation field p, a vector, to describe the 
local orientation of the filaments making up the gel. With every filament one can associate 
a unit vector pointing to one end. The vector p is given by the local average of a large 
number of these unit vectors. Our choice of fiux is {vi,pi), i = x,y. The field p is the 
time derivative of the filament orientational field p. Our choice for fiuxes refiect the lack 
of Galilean invariance in the system. The corresponding generalised forces are (Fj, hi) such 

4 



that the rate of change in the Free energy T per unit surface may be written as 

J d^r[-Vidjafj + h±iPi] = - j d?r[~v,ViTl + htpt] = - j d^r[viFi + hiPi], (1) 

where af^ is the total dissipative stress and we have used the force balance condition djaf^ = 
diU when there are no external forces in the system and Fi = — Vjll. Note, for a pure fluid 
confined between two plates, the velocity Vi in the high friction limit is proportional to the 
spatial derivative of pressure U : Vi oc — Vjll {Darcy's law). In the present problem the 
Darcy's law will be generalised to include forces due to the orientation field. Further, hj_i 
is the field thermodynamically conjugate to pi. Here, " _L " refers to the XY plane. Under 
t — )■ —t, we have Fi — )■ Fi and hi — )■ hi. The rate of change of free energy has, in general 
two parts - reversible J-'rev and irreversible J-'irr parts. The former is the rate of change in 
the internal energy and the latter is linked to the entropy production. The fluxes also, in 
general, can be decomposed into reactive and dissipative parts. Note that the generalised 
forces have well-defined signatures with respect to time reversal. The dissipative fluxes have 
the same signature under time reversal as their conjugate forces, while reactive fluxes have 
opposite signatures under time reversal. In the present problem, both the generalised forces 
are even under time-reversal, and, therefore, the reactive fluxes vanish. Thus, according 



to the Onsager reciprocity theorem (2^, the Onsager matrix will be fully symmetric in the 
viscous limit. Spatial symmetries of the system will dictate the detail form of the Onsager 
matrix. We will see below in Sec. Illllthat our derivation of the 2d constitutive relations from 



the 3d constitutive relations 



19|, |2l|] are fully consistent with the discussions in the Section. 



III. EFFECTIVE 2d DESCRIPTIONS OF THIN CONFINED 3d SAMPLES 

In this Sec. we show how two-dimensional descriptions emerge naturally for thin confined 
three-dimensional {3d) samples with appropriate boundary conditions under averaging over 
the thickness. This averages out the variation along the thin direction and leaves only the 
in-plane variations, leading to an effective 2d description, which may also be derived apriori 
by imposing symmetry conditions between 2d fluxes and forces (see Sees. IVIIII and lIX|) . 

Consider a thin sample of thickness W, confined between two parallel plates, parallel to 
the XY plane. It is assumed to be of infinite extent along the x and y-directions and the 
^-direction becomes confined; W is considered to be much smaller than the lateral in-plane 
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dimensions. In what follows below, we set f ^ = in thin film approximation, which may be 
justified as follows: The incompressibility condition on the 3d velocity field yields 

d,v, + ■ = 0, (2) 

[v_L = {vx,Vy)] yielding (in order of magnitude) Vz ~ x"^-!-' "where L — )• cxd is the lateral 
linear size of the system and derivatives has been replaced by the corresponding inverse 
length scales. Thus in the limit of small W/L, Vz is negligibly small and is set to zero. 
Even though the above approximation rests on the incompressibility of the fluid, we use this 
approximation for the compressible case as well. 

We begin from the 3d constitutive relations and the force balance equations, as reported in 



18[ . The appropriate thermodynamic fluxes are the symmetric part of the stress tensor 



Ref. 

aa/3 and local polarisation field Pa- The conjugate generalised forces are the strain-rate tensor 



Uai3 = {daVj3 + dpVa)l'2. and the local orienting field ha- Here, a, /3 are x, y or z. The rate of 
change of the free energy in 3d = — f d^r[uapcraji = — / dxdy dzluapo-ap + 

VcK] = - f dxdy j^J^i2dz[uci3a^°ji + 'Paha\, where a^°p = cXa^ + cr^^, ci^^ is the antisymmetric 
part of the total stress tensor a^°^. We impose specific boundary conditions (e.g., no-slip) on 
the velocity fields at the confining walls: Va{x,y,z = ±W/2) = 0. This breaks the Galilean 
invariance, as the rest frame of the confining walls becomes the preferred frame of reference, 
in which our equations of motion will be valid. Under partial integrations, the rate of change 
of free energy becomes 

/rW/2 
dxdy dz[-di3aaf3Va + Paha]- (3) 

J-W/2 

In systems with small sizes (biological gels are of the order of micrometers), such as ours, 
momentum conservation is replaced by the force balance condition [l^ 

d^a':; - d^U = 0, (4) 

in the absence of any external forces. Here, a^*^ is the total stress tensor. With the force 
balance condition, Eq. ([3]) reduces to 

/l-W/2 
dxdy dz[-dallva + Paha], (5) 

J-W/2 

suggesting that the pair [va, c^aH) may be treated as a thermodynamic fiux-force in a linear 
response theoretic description. Moreover, treating opposed to flux auto- 

matically breaks the invariance under a Galilean boost: v — )■ v + vq, which is expected in 



the presence of a wall. This is in agreement with our arguments given in SecjlTl Finally, 
averaging over the z-direction and expressing in terms of the mid-plane values of the fields, 
we obtain 

j' = - J dxdy[-djlvi + pih±i], (6) 
which matches with Eq. ([T]). Here, i = x,y. 



The 3d constitutive relations are 18 1 



a = 2r]iu - A/i(^p^ + + ^iPaha (7) 

DPa ha 



+ AiP«A/i - UiUaf3P/3 - VlUPa + ^AP " Vp„ + ^sPaV ■ p, (9) 

Dt 7i 

where a = a^a is the trace of aa/3, d'ai3 = Ca/s — is the traceless symmetric part of 

Similarly u and Uap are trace and traceless symmetric parts of Uai3, rji and 172 are bulk 
and shear viscosities respectively, coefficients z/i, z7i couple with the orientational degrees 
of freedom, ^, ^' and ^ are coupling constants which parametrise the active stress. The 
derivative is the covariant derivative and is defined by -^^^ = Aa+cOapAi^ for any 

vector field Aa with Uai3 = {daVp — dpVa)/'2^ as the vorticity tensor. Terms with coefficients 
and E,B are polar terms which break the nematic symmetry. For a nematic sample 
= = ^B- The antisymmetric part of the stress tensor is given by = ^{pahp —ppha). 
Therefore the total stress becomes 

= 2u6a/s{r]i - ^) + V2{df5Va + daV/s) + {pi - ^)p^h^5ap + ^{Pahp + ppK) 

-A/X(C + i- -^)Sal3p'^ - ^'AfJ,PaPi3 + -iPahp - P^K). (10) 

The force balance condition (jlj) then yields the generalised Stokes equation 

dpa^p = dall = 2{'r]i - ^)9aU + {j^i - ^)da{P'yh^) + m^'^Va 

+ r]2daV ■ V - ^'dp{paPi3)/^^i + ^dp{pahp + ppK) + ]jdii{pahii - ppK). (11) 

We assume there are no mean fiows in the system, i.e., in the steady (without fiuctuations) 
state V = (0,0,0). In addition, we consider three different unperturbed (without fluctu- 
ations) reference states for p with appropriate boundary conditions. In each of the cases 
below, we linearise about the chosen reference states of v and p. We use the Franck free 



energy for the nematic liquid crystals, which, in the equal Pranck's constant limit, is given 

by 

If If rW/2 

^p = 2 J ^^^^^^'^Pl^f ^ 2 J ^^^^ J dzK{daPpf. (12) 

We discard the longitudinal field (equivalently the Lagrange multiplier) as its only role is 
to renormahse the bare coefficients which do not change the results at the scaling level. 
Finally, in some of the cases discussed below we do not consider any anchoring conditions 
on the actin filaments. Apart from its theoretical interests, such no anchoring conditions 
may arise for small-size objects like actin filaments where active anchoring conditions may 
compete against traditional liquid crystal anchroing conditions arising due to energy or 
entropic reasons at the confining surfaces. 

A. System I: 2d Planer polar/nematic order without external forces 

We begin with the reference unperturbed state po = (1,0,0), {vxiVy^Vz) — (0,0,0). 
We show below that the effective 2d descriptions of small fiuctuations around this state is 
a 2d state with planer polar/nematic order without external forces. Since the x-axis is the 
ordering direction, the system should be invariant under y — > —y as nothing distinguishes 
between +y and —y directions. However, there is no x — ?> —x symmetry. Further, since 
Px = I. the fluctuations in p^ are higher order in smallness. Moreover, if the sample is 
nematic, it should also be invariant under p — > —p. This symmetry is however absent 
for a polar sample. We consider both nematic and polar samples here. The effective 2d 
description in this case should be invariant under these symmetries, a fact which we confirm 
below. 

z = w/2 

r r z ^ - " 



z = -w/2 

X 

FIG. 1. A schematic diagram of the chosen initial reference state for System 1. Arrows (parallel 
to the X-axis) indicate the direction of p in the reference state. 
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We impose strong anchoring boundary conditions: Pz = at z = ±W/2 and we impose no 
restriction on py ai z = ±W/2. We first consider the apolar case. Since the fluctuations in p^ 
are second order in smallness we can safely ignore the dynamics of Px ■ Thus the fluctuations 
6p = {0,Py,Pz)- Therefore, we consider the dynamics of Py and Pz only. Further, we impose 
no-slip boundary conditions on the velocity fields i.e. v = at ^ = ±W/2 for all x,y. To 
proceed further, we assume that the dependence of the in-plane flow field on the z coordinate 
is Poiseuille for a given instantaneous in-plane velocity at the mid-plane: We write 

v^ix, y, z) = {z^ - W^/A)e^{x, y), (13) 

which clearly satisfy the no-slip conditions aX z = ±W/2. Here, 6a{x, y) is a vector function 
containing the (x, y) dependences of the 3d velocity fields Va{x, y, z). We assume Py(x, y, z) = 
Py{x,y), i.e. Py has no 2;-dependence to the leading order. Using the boundary condition on 
Pz we can write the simplest form of Pz as 

Pz{x,y,z) = ipzix,y) iz - —j , (14) 

where ipz{x,y) determines the in-plane variation of Pz{x,y,z). In our notations, the fields 
0a{x,y) and ipz{x,y) are proportional to the mid-plane {z = 0) values of the corresponding 
fields, Va{x, y, z = 0) and Pz{x, y,z = 0): 

Vaix,y,z = 0) = -9a{x,y)W^/4: = Vi{x,y), (15) 
Pz{x,y,z = 0) = -^(^x,y)Wy4, (16) 
Py{x,y,z = 0) = py{x,y), (17) 

since Py is assumed to be independent. Here i = a = x or y and Vi{x,y) is a 2d velocity 
field. With the choice of ^-dependence of Py and Pz, the Franck free energy J-p in 3d becomes 

^P = \j dxdyK[^pl{x,y) + W{diPyf]. (18) 

Clearly, from the expression (fT8|) above, Pz is massive and does not survive in the long time 
limit, i.e., pz is not a slow variable. This is due to that fact pz is held fixed at the walls at 
z = ±W/2. In contrast, Py is a massless field and its fluctuations are long lived; Py remains a 
slow variable in the problem. Further, Eq. flTH]) allows us to define the 2d effective conjugate 
field h±i as 

h^. = (19) 

OPi 
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Now we average the force balance equation Eq. (fTTj) over the thickness W to obtain an 
equation as a function of x and y only and are written in terms of 2d vectors. After 
averaging, in the lubrication approximation 22| we get 

T77 V2y^v^dz = --^v,{x,y). (20) 



After performing 2;-averaging over other terms, we find for the y-th component of the velocity 

dyli = -^Vy - ^'d^{p^Py)Afi + ^^^^ 9^ ± y ) (2l) 

-Ddyll - ^2dxiPxPy)AfJ' + i^2dxiPxh±y), (22) 



u 



y 

where constants D = the inverse friction coefficient, ^2 = ^'D and z/2 = ^f^-D. Similarly, 
for the x-component we get after averaging over z 

= -yy^'^x - Cdy{pxPy)Afi + -^^^dyipxh^y) (23) 
^Vx = -DdxU - ^ody{PxPy)Afi + vody{pxh^y) (24) 

where Vx and Vy are components of the 2d vector t>j; = i'D and vq = ^f^-D. 

Having obtained the effective 2d Eqs. of motion for Vx and Vy above we now proceed to 
obtain the same for p. Linearising about the reference state chosen and neglecting the terms 
higher order in smallness Eq. reduces to, 

^ = — - ^^ {dyVx)Px - ^ {9xVy)Px + U^xPy (25) 

ut Z Z 

Averaging over the thickness W , we get 

dpy _hy '^l + 1 1^1 - 1 



dt Jo 3 
where "Jo = liW and Px = I- 

From Eq. ([22]) and Eq. we get 



idyVx)Px idxVy)Px + U^xPy , (26) 



dxVy = -Ddxdyll - ^2dl{pxPy)Afi + V2dl{pxhy), (27) 
dyVx = -Ddxdyll - ^ody{pxPy)Aij, + uodyipxhy). (28) 

Putting Eq. (^3) and Eq. (^E^ in Eq. (^E^ we get 

dpy h±y 2pxWvl^ 2pxWvl^ '^pIWvq ^ 

-Qf = ^ ^^<9j^(Px/i±y) ^j^d,{pxh±y) + 3^ ^oAfidyipxPy) 

+ ^^^^2Af,dl{pxPy) + (z/0 + iy2)^^dxdyU + UdxPy. (29) 
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Clearly when A/i = 0, Eqs. (122|24|29p do not immediately conform to the (symmetric) 
structure as stipulated by the Onsager Reciprocity Theorem. In order to achieve that we 
exploit the rescaling freedom of the fields that comes due to the (multiplicative) arbitrariness 
in the definition of 2d py. Let us scale Py by A where A is any real number. Then py — Xpy 
and hj_y — Xh±y. Scale factor A is to be chosen such that for AyU = 0, a fully symmetric 
structure of the coupled equations for Vx, Vy and Py follow. The scaled equations are 

Vy = -Ddyll - A^2<9x(PxPj/A/i) + \i^2dxiPxh±y), (30) 

Vx = -DdxU - X^odyipxPyAfi) + \uody{pxh±y), (31) 
dpy h^y 2pxWvl^ 2pxWvl^ 

^ = — ^dy{pxh,y) - -^dAPxh.y) 

+ ^-^^,A,.dl{PxPy) + ^-^i.^l^dliPxPy) + (.0 + ..)^5.9,n. (32) 

We choose \ = {^y/'^. Thus Eq. (El becomes 



h^y 2W 



yldl + yldl)h^y + (^0 + z/2)(^)^/'a,a,n 



+ |^A/i(z/oeo5j + z/265,2)p, + UdxPy. (33) 

Now we rescale the coefficients (^)^''^(z^2, ^'o, ^2-, ■Co) (^^2, J^o, ^2-, ^o)- Hence, Eqs. ( l30]l3Tll 
and ( l33l) become 

Vx = DFx + vodyh^y + ^oAfidyPy, (34) 

Vy = DFy + U2dxh±y + ^2^fidxPy, (35) 

-^^oAfld^Py - ^^2A^idlpy + ^A^xPy (36) 

It is evident that effective 2d Eqs. fl34f35p and (!36|l have a symmetric structure for AyU = 0, 
as expected. Further, they follow the same symmetries as discussed at the beginning of the 
present section. An independent direct derivation of these 2d equations, based on symmetry 
arguments, is provided in Sec. IVIIII 



B. System II: Polar order normal to the plane without external forces 



We take the reference unperturbed state as po=(0, 0, 1). Thus the ordering is normal to 
the surface of the plane. Therefore, one has in-plane rotational symmetry in the system. 
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We will see below the small fluctuations around this state allows for a 2d descriptions with 
polar symmetry without any in-plane polar order (the macroscopic polar order lies normal 
to the plane) and without any external forces. As before, we set Vz = 0. Further, in the 
fluctuating state (p^) = = (py) and (v^) = = (vy) (since there is no external force); thus 
Pi and Vi are of the same order of magnitude. 



z=w/2 

I ! I t t 1 
t I t t t ! 
t t I t t t 

z=-w/2 



FIG. 2. A schematic diagram of the chosen initial reference state for System II. Arrows (parallel 
to the z-axis) indicate the direction of p in the reference state. 



Since = 1 in the unperturbed state and we are looking for small fluctuations around 
it, fluctuations in pz are second order in smallness. Hence, the lowest order polarisation 
fluctuation 6p= {Px,Py,0)- Henceforth the dynamics of Pz is ignored. We consider the 
dynamics of p^ and Py as both are slow modes. We assume p^ and Py to be independent of z 
to the leading order and impose no anchoring conditions on them. Next we impose no-slip 
boundary condition on the velocity flelds: Va{x, y,z = ±W/2) = 0. We impose the thin fllm 
approximation on the generalised Stokes Eq. ffTTj) and average over z. We flnd, as before, 
W /-l^/2^2V^fa dz = -^Vi{x,y). In addition, 

^ l'W/2 i-W/2 

TTz f u dz = - I dpVpdz — )■ subleading. (37) 

All the other terms, upon 2;-averaging, turn out to be second order in smallness. Therefore, 
new terms are to be added in the ?>d onsager relations. The next order terms are polar terms 
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and we have, 



cr = 2r]iu - A/i(^p^ + + ^iPaha + ■ h - .C2A/iV ■ p, 



(38) 

(39) 
(40) 



71 

We henceforth set z/i = z/i = ^ = ^ = ^' = 0, since they do not contribute anything at 
the hnear order. The modified force balance condition yields: 



1^2 



^2 



a„n = 772V Vo, + (z/2 - -7)'9aV ■ h + —dfi{djifi + df^h 



-(6 - -|)A/XaaV ■ p - -^AfxdisidaPfS + dfsPa)- 



(41) 



1 



Eq. ( HTl) has a term Y'^l^i^ahp + S/j/iq.) which is equivalent to ^dj{dihj + djhi) as 
and = 0. Here i,j = x,y only. In contrast pa has a term ^dpUap which comes out to be 
Y^d'^Va + 2djUaj)- Hence, to keep the force balance equation and p'a on the same footing we 
modify the cross coupling term in the pa Eq. to U2djUaj, (j = x,y). After the z averaging 
Eq. pT]) . the equations of motion for the 2d velocity fields fj are obtained: 



Vi = -DViU + ^dj{djhi + dihj) + V2odiV ■ h - ^AjjdjidjPi + dipj) - ^2odidjPj, (42) 



where = ^20 = i>2D, U20 = {1^2 - f )I^, 60 = (6 - %D and 60 = ^2^- 

After z averaging the third of Eqs. (HOl) resulting 2d equation of motion of pi at the linear 
level is 

Pi = V —djUij + -[1J2- -rjOiUi, (43) 

where 70 = 71 TV. 

In Eq. (HHj) we substitute for fj by using Eq. to obtain a constitutive relation for pi. 



Pi = \ —^±h±i + 

7o 6 



I/2l^20 , 2l>2Z/20 , 2z/20, Z>2 x , ^ ^20 ^^2 . 



V,;Vi V 



^2-D„9 



Z>2/^ 2D,_ i)2. 



ViVi ■ F - ^e2oA/iVipi 
o 



'^2'C20 _|_ 2 z/2 ^20 



A//ViV,V± ■ p - ^(z?2 - ^)(6o + 6o)A/iV,ViV± ■ p. 



(44) 



2 3 

Note that Eqs. (H2l) and (H4|) do not really conform to the Onsager symmetry: In order to 
make the system in agreement with the Onsager reciprocity theorem, we make the rescaling 
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Pi — )■ Pi\/2/\/3. The final equations for Vi and Pi take the form: 

= -DV.n + ^mK. + d.h,) + ^5.v^ ■ h - ^A^^^,i^,p, + d,p,) - ^v,v, 

Pi 



+ ^dMF^ + d.F,) + ^V.V,F, - [Vi,. + 3ViV.V. ■ p] 

- ^A;.V.Vl V. ■ p - ^(6o + 6o)A/.V.Vl V. ■ p. 

In Eq. ( H3|) one could have added another fiow-orientation coupling term proportional to 
AfiVi which would have originated from a term of the form Afid'^Vi in the 3d Eq. for pa- 
Such a term, being proportional to A/i, vanishes in the equilibrium limit and hence keeps 
the Onsager symmetry unchanged. When A/i 7^ this term generates a leading order active 
term AfidiU in Eq. (jH]), in addition to several subleading terms. Such terms however, being 
products of AyU and diU, beyond the scope of linear response regime. 

By using the symmetry arguments given at the beginning of the section we can arrive at 
a similar set of equations directly, instead of averaging the 3d equations over the thickness. 
What results are equations identical to Eqs. fH^ . 

Vi = DFi + P2^^{y± ■ hx) + i/sVi/i, + e2A/iVipi + ^sA/i^iVx ■ p, (46) 



Pi 



7 D -L J J. J J. J 



h±j + U2^^iV± ■ F) 



+^3ViF, + i^^i±Me3A/.Via,(Vx ■ p) + ^6A/xVi9.(Vx ■ p) + ^A/iVip.(47) 

Now comparing the above equation with the one obtained from Eqs. ( 145|) we get the 

following relationship among the various coefficients in the two different cases: ^ = 1^3 = 
T, — I fid \ r, \ <: — €20 f _ \pi ( C20 



C. System III: 2d planer polar order with external forces 

In this active gel film the Id symmetry present in the sample is y — —y. Since p^. = 1, 
there is no invariance under x — )■ —x and the system being polar it is also not invariant 
under p — t- —p. Furthermore, Vi is no longer a gradient due to the presence of external 
forces. The effective 2d equations of motion of this sample should be invariant under these 
symmetries. We write down an effective Id description by averaging over the 3d sample. In 
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the above two cases we considered above the chosen reference states were uniform. However, 
we now choose a non-uniform reference state. The reference state is given by, 

= 0, = ±1, py = Oatz = ±W/2 
= 1, = 0, Py = at z = 
Py = everywhere in the bulk (48) 

We parametrise the initial state po by 

p^ = cos 9{z), p^ = sin6'(2:), Py = 

Before discussing fluctuations about this reference state in details, let us point out the 
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FIG. 3. A schematic diagram of the chosen initial reference state for System III. Arrows indicate 
the direction of p in the reference state. 

difference between System III and Systems I and II. In the latter two, the transformation 
Px — )■ —prc either makes p — ?■ — p or p — )■ p in the reference state: In particular, this 
transformation for System I makes p — — p (since Pz = Py = and for System II makes 
p — 7- p (since Px = Py = 0)- Further, since microscopically the molecules are nematic, 
J'3d{Px) = J^3d{~Px)- In contrast, the transformation p^ — )■ —px, when applied on System 
III, does not make p — )■ —p. See Fig. HJand Fig. |5] showing System III with p^ — )• —p^ and 
p —7- — p, respectively. Clearly, Fig.HJand Fig.[5]are not the same. Nematic symmetry implies 
J^3d{p) = -^3d(— p)- For Systems I and II, this immediately implies J-'sdiPx) = J^3d{—Px)- 
However, for System HI, J^sdiPx) 7^ J^3d{—Px)- Since effective 2d descriptions are constructed 
in terms of the mid-plane values of the 3d vectors, we expect J^px 7^ J^-p^ , indicating that 
the effective 2d descriptions for System HI is expected to show up polar symmetry, which 
we confirm below. The Franck free energy in the equal Franck's constant limit becomes: 
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FIG. 4. System III with 
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FIG. 5. System III with p 



dxdy 



K 



W/2 



W/2 ^20 

dxdy I 9(z) 

w/2 



dz'^' 



for = cos6'(z), = sm6{z) and = 0. The Euler Lagrange equation reads 

d^e 



dz"^ 







6 = Az + B. 



(49) 



(50) 



which thus defines the reference 



Imposing the boundary condition we obtain 6{z) 
state. 

We impose no shp-boundary conditions on the velocity fields. There is no restriction on py 
(only Py is the slowest mode) at the boundaries. Then z-averaging of the generalised Stokes 
Eq. f|TT]) together with the assumed 2;-dependence of the velocity field given by Eq. f|T3l) we 
obtain 



-Ddyli + l'2h^y + C.2^f^Py, 

-DdxU - uodyhxy + ^o^f^dyPy, 



(51) 
(52) 



where = D 



and 6 = ^, .0 = and eo = fi^. 
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Linearising Eq. (Q for the y-component about (p) = cos6{z)ex + sm6{z)ez and averaging 
over the z components we get 

OtPy = h 8 Vy + 8 :^dyV^, (53) 

where 70 = 71 T^. Substituting for and in the above equation and rescahng = 
-^a,n and Fy = -'^dyU, we get 

= D^F^ + uodyh±y + ^oA^dyPy, (54) 

Vy = DyFy + J_y + ^ ^/^P^ , (55) 

-uodyF^ - ^^oA^dlpy. (56) 



where 7 = 'yiW,Dx = Dy = D. In Eqs. and ( 1551) although we have = Dy, in 
general 7^ -D^. See Sec. IIXI for an alternative derivations, based on the 2d dymmetries, 
where such anisotropic coefficients appear naturally. Finally, one may generically add polar 
self-advecting terms p ■ Vp and pV ■ p in the equation for py in Eq. (!56|) . Upon linearisation, 
we find 

1 z/2 
= [- + ;^ - -^9l]h^y + u^Fy - UdxPy 

-uodyF^ - ^^oAfid^Py, (57) 



where is a coefficient of either sign. Equations (jSH - [57]) formally constitute effective 2d 
description for System III. 



D. Dynamics of actin filaments in a motility assay 

Motility assays are standard experimental set ups to measure properties of molecular 
motors such as myocins. In such an assay one observes in the light microscope the motion 
of isolated actin filaments being propelled by the myosin heads which are immobilised on a 
glass surface. Actin filaments are observed to slide over a layer of motor proteins (myosin) 
bound to the surface (glass plate). Actins bind with and unbind from myosin heads stochas- 
tically, resulting into velocity (or momentum) being imparted to the actin filaments. The 
actin velocity has a mean and fiuctuating parts. It has been observed that actin velocity 
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increases with ATP concentrations 25|. In recent controlled experiments with many fil 



aments at high densities, complemented by cellular automata based simulations |26|, the 
authors demonstrated the spontaneous emergence of order and collective motion in the sys- 
tem. Further, they find large density fiuctuations and propagating modes in the system. In 
this Sec. we discuss a simple coarse-grained description of the coupled dynamics of local 
actin velocity and polarisation degrees of freedom in a motility assay. We continue to use a 
one-fiuid model, and thus do not distinguish between the solvent velocity and the velocity 
of the filaments (with the velocity field in our model equations corresponding to the centre- 
of-mass velocity in a two compoment system). This, although clearly an approximation, it 
suffices for our purpose to show the generic polar nature for our system and its connection 
(in the sense of same symmetries) with our System III discussed before. We assume that 
the actin filaments are preferentially oriented along the a;-direction, i.e., (p) = 1. Assum- 
ing incompressibility, we begin with the generalised Stokes Eq. ( ITTj) . However, in an actin 
motility assay, the velocity at one (or both) the walls is not zero, due to the forces imparted 
by the motors bound at the walls. Thus, instead of no-slip boundary conditions there is a 
finite slip velocity at the walls. Without any loss of generality, we assume that the surface 
forces exist only at the wall at z = W/2. In the presence of the external forces, the force 
balance equation (jl]) is generalised to 

d^a':^ - d^U = /r , (58) 

where is the total external force. These external surface forces, on symmetry grounds, 
should depend on either ha or pa, and hence have two contributions of the form (i) haS{z — 
W/2) and paAfi6{z— W/2) in the hydrodynamic limit. The 5-function ensures that the forces 
exist only on the top surface at z = W/2. The first contribution survives in equilibrium, 
while the second one is essentially a nonequilibrium contribution. 

Together with these additional forces, the generalised Stokes equations for Vx and Vy, that 
follows from the force balance equation flSSl) . for an actin motility assay become 

r^VV. = -Dd^U + C.'dj{pxPj) - ^d^ip^h^ + p^h^) - ^^dpip^hp - pph^) 

+ Dohx6{z - W/2) + CoPxA/i5(^ - W/2), (59) 

r]V'^Vy = -Ddyll + i'd^{pyp^) - ^d^{pyh^ + p^hfj) - -dpipyhfi -pphy) 

+ uohy6{z - W/2) + CoPyA/x5(z - W/2), (60) 
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where, coefficients uq and Co couple the velocity fields with the surface force terms. Averaging 
Eq. (!60|) over the thickness and in the lubrication approximation, we obtain for the x and 
y-components of the 2d in-plane velocity field Vi 



U 



y 



-bdjl + i>o/i_L^ + i>xdy{pxh±y) + ioAnp±^ + i^dy{p^PyAij,), (61) 
-Ddyll + uoh±y + Uydx{pxh±x) + ioA^pj_y + iydxiPxPyAfi), (62) 



with coefficients ^o, ix, iy, i^o, ^x and Uy which are related to the coefficients in Eqs. ( l60l) : D 
is an inverse friction coefficient. To the leading order in gradients equations (16111621) have 
exactly the same form as those of Symmetry III in Sec.[Tll Clearly, velocities are proportional 
to A/i and the mean velocity (vx) increases linearly with A/i, or, equivalently with the ATP 
concentration. Since in this system {p) = poix = ix to the leading order, we ignore h±x- 
Therefore, 

Vx = -DdxU + i>xdy{pxh±y) + ^oAfipx + ^xdyipxPyAfi), 

Vy = -Ddyll + i'oh±y + ^^A^p^y + ^ (Px'Pj/ A /i ) . (63) 

Here, ^qA^p^ gives a constant contribution to Vx- This term is responsible for the steady 
motion of the actin filaments as observed in actin motility assays. By comoving we can get 
rid of that. Therefore to the leading order, 

Vx = DF^ + Uxdyipxhxy) + ixdyiPxPyAfi), 

Vy = DFy + z>o/i_Ly + ioAfipy, (64) 

where we have used Fj = — Vjll. Equations (l64l) are identical to the constitutive relations 
set up for Symmetry III above. 

The 3d bulk equation for the polar order parameter in an one component description is 



of the form 



2l| 



+ V ■ Vpa + UJapPp + ^AP ■ Vp + ^spV ■P = ^ + AVn + VlUadPli, (65) 



where Ho is the pressure 2l| and A is a coupling coefficient. To proceed further, we linearise 
Eq. (|65l) about a given direction (say x-direction) for macroscopic orientation, we obtain an 
equation for py of the form of Eq. ( 1571) . This shows the close connections between our System 
III and the coarse-grained dynamics in an actin motility assay. More complete analysis of 
actin dynamics in a motility assay would require a two-fiuid approach and proper boundary 



conditions on the free surface which we defer for a future work 

19 



27|. 



IV. STOCHASTIC DYNAMICS AND TIME REVERSAL INVARIANCE 



The invariance of the statistical steady states under the time reversal operation charac- 
terises systems in thermal equilibrium. An important consequence of this is the Fluctuation 
Dissipation Theorem (FDT) [l| which connects correlation functions and corresponding sus- 
ceptibilities in equilibrium. For systems outside thermal equilibrium, there are no FDT. In 
order to calculate correlation functions from the Eqs. of motion set up above, we add noises 
which represent the inherent stochasticity of the models arising from the fast degrees of free- 
dom which are ignored in the dynamical descriptions of the models. Stochastically driven 
coarse-grained models have a long history in statistical mechanics, starting from equilibrium 
critical dynamics [3] and later on in nonequilibrium phenomena. Notable examples of the 



latter include surface growth phenomena 



driven diffusive systems 



29| . diffusion medi- 



ated reactions [30] etc. In such stochastically driven models added noises are assumed to 
be zero-mean Gaussian distributed with specified variances. For systems in equilibrium, the 
variances are fixed and are linked to the dissipative kinetic coefficients of the model. This 
arises as a consequence of the FDT. In systems of biological interests, the noises are of both 
thermal and non-thermal origin. In this article, we confine our discussions to the thermal 
noises whose variances can be fixed by using the FDT for A/i = 0. We work out explicitly 



for System I below. Introducing thermodynamic forces 
Onsager relations can be written in a matrix notation as 



—dril and F„ 



—dyll the 



If \ 



Li 



\PyJ 



( 



Li 




D 



\ 



1 2 



(66) 



D 

U iV2qx 
\^ -iv^qy -iv2qx ^ + ^ + ^ ) 
Hermiticity of L/ follows from the fact that all the forces F^,, Fy and h^y have the same 
property (even) under time reversal [l^. The time reversal property of the force can be 
found out from its definition. For a dissipative flux, its signature under time reversal is 
same as that of its corresponding force. Positivity of dissipation ensures that the dissipative 
coefficient relating a fiux-force pair is strictly positive. When A/x = the noise correlation 
matrix = ^k^TLi (see below). Such a choice ensures the FDT. When A/x = all the terms 
in the and Vy equations are even under t — )■ —t. Since the active terms are proportional 
to they are odd under t — t- —t. Thus AyU, a nonequilibrium drive, breaks the symmetry 
under the Onsager reciprocity principle. In order to show that our choice of the noise 
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variance matrix L indeed satisfies the FDT, we calculate the susceptibility and correlation 
matrices for the dynamical variables in the presence of external sources. For System /, the 
relevant dynamical variables are density p and orientation field Py. Their coupled dynamics 
may be written in a matrix notation as 

/ 



p \ f -1 / p \ , f / no 



. -^XoM ' I, (67) 

\Py ) \Py J \ h±yO 

where Ho and hyo are externally imposed forces. Matrix xo is the static susceptibility matrix 
and in the Fourier space is given by 



Xo ^ = ' : ■ (68) 




Matrix L in the Fourier space is given by 



L 



(69) 



The dynamic susceptibility matrix x is given by x = Xiji^) = [i-^^ + Lxo] ^L. Here / is the 
unit matrix. For a single-field model, the FDT yields a relation between the susceptibility 



x(a;) and correlation C{uj): x"(a;) = |a;C(a;), P — 1/KbT where x"(^) is the imaginary 
part of x(^)- Foi" ^ multi-component theory, such a relation generalises to a matrix equation 
and holds element-by-element: x'lji!^) = 2^Cij{u), where i,j refer to different fields (here 
p, Py), Xij is the imaginary part of the susceptibility matrix element Xiji^)- To calculate 
the correlation function matrix C = Cij{uj), we add thermal noises fp = —iqxfx — 'iQyfy aiid 
fp, respectively, in the Eqs. for p and Py. Here, fx and fy are the thermal noises for Vx and 
Vy. Noting that the variances of the noises fx, fy, fp are given by 2KbTLj, we obtain the 
variances of the noises fp, fp which is just 2KbTL. Matrix C_ is then given by 

C{uj) = 1^ ^ j (p Py)* = [iuji + Lxo\-'2KBTL[-iujI + XoL]''- (70) 

Further, from the definition of Cij{uj) we have Cij{uj) — Cji{uj) yielding Xiji^) — ~Xji{~^)- 
We, therefore, conclude that Xij is the anti-hermitian part of the full susceptibility matrix 
X- One then obtains 

X" = [tcol + Lxo]-'ujL[-iujI + xoL]-' = ^ujC{uj), (71) 

estabhshing the FDT for the present case. This shows that our choice of the noise variance 
matrix is in agreement with the FDT. 
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V. LINEAR STABILITY ANALYSES AND CORRELATION FUNCTIONS 



In this Sec. we set up the Eqs. of motion for the relevant slow variables, starting from the 
2d constitutive relations which we established above for the different cases. We separately 
consider the compressible and incompressible cases of the dynamics. 



A. Dynamics for System I 

We now set up the equations of motion of the slow variables for System I, for both 
nematic and polar orders. The slow variables are Py and mass density p when the system is 
compressible; in the incompressible limit p, being a constant, drops out and Py remains the 
only slow variable in the problem. When Ap — 0, the equations of motion are consistent 
with the Onsager reciprocity principle. The continuity equation for density p reads 

|e ^ _V . (PV) = -mpV,) + dy(pVy)]. (72) 

Now linearising about the mean density (p) = po = 1 we get 

^ = DVlU - (uo + V2)d.dyh^y - Ap(6 + io)da.dyPy. (73) 

The free energy J^p for a polar sample, in the equal Franck's constant limit in 2d {ki — 
K2 — k), is given by 



J 7 



p = y dxdy{~p^\pi + ap • V_lP + \^XpP^\ = J dxdyfp. (74) 

Parameters a and Xp ^-re taken to be constants; a couples density fluctuations p with 
polarisation p, Xp for a sample with p = is the inverse compressibility. For a nematic 
sample J-'p must be invariant under p — p and hence a = for a nematic sample; for a 
polar sample a 0. 

From thermodynamic considerations we now use (linearising about a mean density po), 

n = pj^ - fp^ iXP - aV± • p), and h±y = ^ {i^^lPy + aV±p) (75) 

Therefore, Py in the Fourier space is written as. 



Py = [— ^ ^ H^o + i^2)aqxqy + —^o^pQy + —^2Apq^ - iUQx\Py 



+ + + - (^0 + V2)qxqyX\P- (76) 
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Equations (173|) and (176ll have been considered in Ref. [Sj in a discussion on active nematics. 

The equations of motion for py and p can be cast in a matrix form whose eigenfrequency 
in the hydrodynamic hmit (small q) is given by 

2 ,/„J^„A,,«2 



Afe,g,) = ^[-Dql -■^ + ^^^V^ + ^-^^^^ ± {{-Dql + 



y 



2' 7 D D 7 D 

'-y - 4(eoi^o + ^01^2 + ^2^0 + i2y2)^MY'\ (77) 



D 

We can now compare eigenfrequencies ( !77|) with those in Ref. [8|]. Although our Eqs. for 



System I are invariant under the same set of symmetries as those in Ref. 



there are 



some differences in details. For example, the control parameter a for creating instability in 
Ref. js] appears in the cross coupling term in their density Eq.; in contrast, in our model 
A/i plays similar role and appears in the orientation equation direction, in addition to the 
density equation where it appears as a cross-coupling coefficient. As a result, our model 
has more complicated behaviour (in regard to the presence of propagating modes, damping 
or instabilities) as a function of A/x, although the general scaling behaviour is identical to 
those in Ref. j8|. One may further calculate the two eigenvalues separately as functions of 
(Q'x, % = 0) and {qx = 0, qy) and analyse their stability: 

= = "-2^ - ,U,., -Dix, , (78) 



1 



KQx = 0, qy) = ^^1-^^^% - l^^oly + l^oCo^my - D 7xg; 



±{{DKql + n^ulql - 7z/o^oA/igJ + D'^^x^lf 
-AD-f{-DWq^y - Da^vl^ql + D^nxqt, + Dn^vlxql 
-D^uoio^y^Xqt)Y'^- (79) 

Equations flTHl) and fl79l) shows that depending upon the sign of the coefficients one or 
both the eigenvalue(s) may become unstable. For '^^^^^^ > and > |^| one of 

the eigenvalues \{qxi<ly = 0) changes sign and leads to low wavevector instability of the 
corresponding eigenvector. The other eigenvector remains stable. Therefore, physically only 
a specific linear combination of Py and p, determined by the unstable eigenvector, is unstable. 
However, there are no instabilities at q = as both the eigenvalues vanish. Further, for 
large enough q, both eigenvalues are negative. In case of equilibrium {Ap = 0) both the 
eigenvalues are stable. 
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Having considered the linear instabilities, we now investigate the opposite situation when 
the system is linearly stable and admits nonequilibrium steady states. This happens when 
the condition of linear instability is not satisfied, i.e., when the active stress is either tensile 
(AyU > in our notation) or it is still contractile, but |A/i| < |AyUc|- The correlation function 
Sq = {p{(l)p{—(l)) in the hydrodynamic limit can be calculated adding noise terms in the 
equation for p and Py. The Eqs. (!73|) and (!76l) take the form 

dtp = xDV'^p - (^0 + ^2)^P'dxdyPy + ViO, 

dtpy = T{q)q'^py - (z/Q + V2)qxqyp + fp- (80) 

where Q is a Gaussian-distributed thermal noise for the velocity degrees of freedom with zero 
mean and a variance (0(x, t)Cj(0, 0)) = 2D(5(x)5(t) andr(q)g2 = + i^a^g^ ^ i^2^g2 

We have considered the case when ^a = ^- From these equations we obtain the equal-time 
correlation function (or structure factor) Sq = {piq)pi—q)) = f ^{piQ,^)pi—<l,—^))- Struc- 
ture factor Sq may show different behaviour depending upon whether or not propagating 
modes are present. For the choice of parameters when there are no propagating modes, we 
find after a straight forward algebra 

Sq ~ (81) 

suggesting density fiuctuations diverge as for a 2d system with a linear size L j^. In 
an equivalent equilibrium system this result would mean the compressibility x diverging as 
in the g — )■ 0, or as in real space, where L is the linear system size. Since ~ A^, 



where N is the total number of particles, density fiuctuations diverge as (see Ref. |8l. l20| ) . 
However, when there are propagating modes (or when 7^ and the system is polar), such 
strong fiuctuations get cut off. 



1. Incompressible limit 

In typical biologically relevant situations, fiow velocities are of very small magnitudes 
(typically much smaller than the sound speed). As a result, the system behaves as an 
incompressible system and we set density p = 1. Incompressibility is enforced by V ■ v = 
dxVx + dyVy = 0. Using this we can express pressure in terms of the other remaining fields 
in the problem and eliminate it from the dynamics of the problem. We have 

H = -^^[(t'o + i'2)dxdyh^y + (^0 + ^2)^P'd:,dyPy]. (82) 
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Substituting for 11 from equation (182|) . we therefore obtain for py in the Fourier space (in 
the polar coordinates Qx = q± cos 9 and qy = q± sin 9) 

Py = - —q (z/Q sm 9 - U2 cos 9) py + -—q (z/242 cos 9 + uq^q sm 9)py 

'J JJ L) 

^^^K^'o^ + ^^2^0) sin^ 29py - iUq cos 9py. (83) 



4D 



Eq. ( l83l) is similar to the one derived in Ref. SjJ for fixed height and constant density. 



Writing Py{q,t) ~ exp[A(g, 6')t], we obtain for the eigenvalue A 

A(g, 9) = - ^q\uo sin^ 9-^2 cos^ 9)^ + ^g'(z/26 cos^ ^ + z/o^o sin^ 



(^^06 + ^^2^0) sin^ 2^ - i^q cos ^. (84) 



AD 

Clearly for Afi = there are no instabilities which is expected and, with = 0, equilibrium 
decay results are recovered, such that A(g, 9, A/i = 0) < 0. For non-zero A/i, two distinct 
cases arise: (i) When X{q,9) < 0, i.e., the system is stable, and (ii) when X{q,9) > 0, i.e., 
the system is unstable. In the stable case, the role of A/i is to enhance the decay rate of 
fluctuations. Further, there are propagating modes of speed C,a(1cos9. Moreover, even when 
there are instabilities, the real part of A(g, 9) initially rises from zero, reaches a maximum 
at g = qmax and then decreases to eventually become negative. We find generally 

D , A/i 

27(z/o sin^ 9 — U2 cos^ 9Y 2k[uq sin^ 9 — 1^2 cos^ 9y 

)i/2. (85) 



z/2^2 cos^ 9 + z/q^o sin"^ 9 - ^(z/06 + ^^2^0) sin^ 29 



A schematic plot of the real part of the unstable eigenvalue as a function of wavevector is 
shown in Fig. |6l It reaches a maximum and then comes down to zero to become negative 
eventually for high wavectors. Since the instability peaks at wavenumber qmax, in the cor- 
responding real space picture there will be a pattern periodic at a scale ~ 1/qmax- Since 
Qmax depends upon 9, the generated pattern is clearly anisotropic. A schematic plot of qmax 
versus 9 is shown in Fig. [3 We clearly see (i) qmax depends strongly on 9, an illustration of 
the ensuing anisotropic pattern and (ii) for larger A/i(A/i// > A/i/) qmax is larger, implying 
that as A/i rises patterns become denser in the real space. 
We analyse for a few different values of 9. At 9 = n/A 

Py = ^ - -^('^o - '^^YPy + ^^'^'^26 + '^o^o - ^^06 - ^^2^0)^3; - -^qUPy (86) 
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FIG. 6. A schematic plot of the real part of unstable eigenvalue A versus wavector q. The vertical 
broken line gives a measure of Qmax ■ 
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2.5 



FIG. 7. A schematic plot of qmax versus 9 for some values of the parameters in System I for two 
values of A/i : A/.i// > A^/ (see text). 



In general however if 



l-g' - ^^{^2^2 + ^^o^o) + ^{^oi2 + i^2io)q^ + ^(i^o - ^2f\ < 0, (87) 



then there would be long wavelength instabilities, with q^ax = l[^ ~ ^{^2^2 + t'o^o) + 
^(^^06 + 1^2^0)1/4^(1^0 - T^2?\. Next, consider ^ = 0. Then, 



Py 



-Py T^Py 



-Py - iUqi-Py 



(88) 



7 - " D D 
Thus there would be long wavelength instability if ^^^'"^ > and |52^/^| > ^. Threshold 



of instability is given by ^^-^ 



0. At 6* = f also we get the same result as in 6* = 0. In 
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general, since A depends explicitly on the polar angle 6, it is entirely possible to have stable 
and unstable regions coexisting together for different values of 6. For example, for angles 
determined by the inequality (with ^o^/^ > and ^2^/^ > 0) 



+ — (i^o6 + ^^26) A/i sm^ 29 < 



7 



4L) 



D 



V2i2 cos^ Q + vq^q sin^ 



(89) 



there will be instabilities, where as for angles outside this zone, it will be stable. On the 
other hand, for sign reversal of ^o^/^ and ^i^ii, 6 satisfying Eq. fl89|) corresponds to stable 
region and 9 outside this domain corresponds to instability. Thus, depending upon the 
parameter values, certain anisotropic (i.e., ^-dependent) patterns will emerge in which py 
will have large magnitudes for certain values of 6', whereas Py will decay to zero for other 
values of 9. In each of these cases, replacement of the inequality sign by the equality sign in 
the relation fl89|) yields a critical activity AyUc above which the instability sets in. A contour 
plot of A versus q and 9 is shown in Fig. [8] clearly depicting negative (stable) and positive 
(unstable) regions. 

In terms of the relations between the different coefficients as shown in Sec. IIII At we can 
rewrite the instability condition fl87|) . which shows that there are different critical thicknesses 
along different directions of the Xy-plane. In particular, we find for = the critical 
thickness Wco is given by 



W, 



cO 



1/2 



_7i(z/i - l)e'A/iJ 
and for 6* = 7r/2, the critical thickness Wctt/2 is 



(90) 



W. 



c-k/2 



nl/2 



7i(z/i + l)eA/i 



(91) 



Since Wco 7^ Wcn/2, instabilities set in at different thickness in different regions of the XY- 
plane for a given A/i [see Fig. (|8])]. 

The correlation function in the incompressible case has to be (p^(r,t)) as Py is the only 
relevant flux in the system. We add a noise ipp to Eq. ( l83i) whose correlation is chosen so as 
to ensure FDT for AyU = 0. 



{ijp{q±,u)i:p{-q±, -u)) = 2kBT 



1 g^(z/o sin^ 6^ — cos^ 6')^"' 
7 ^ D 



(92) 



Correlation (p^(r,t)) reads (p^(r, t)) = f^72^(|Piy(q±, i^)P)- Now from this, using the 
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FIG. 8. (Colour online) A schematic 3d contour plot of the eigenvalue A as a function of wavevector 
q and polar angle 6 for some specific parameter values for System I. Both positive and negative 
values of A may be seen (see text). 



value of (l^/'pp), we get 



{pl{v,t)) = ksT j 



X 



7 



+ 



q±dq± dO 
2tx 27r 

4 



- + ^(z/osin^2^-z/2cos2^)2 
7 D 



D 



vq sin^ 6 — U2 cos^ 



Au 



cos"^ 6^ + i/Q^o sin"^ ^)gi 



(93) 



If - if (^^26008^^ + z/oeosm^^)gi + ^(z/q^ + z/2eo)gi sin^ 2^] > for all 9, the sys- 
tern is stable and {py{r,t)) has a logarithmic infrared divergence. In contrast when — 
^('^26cos^^ + z/o^osin^6')gi + ^(z^o6 + ^^2^0)?! sin^ 261] = 0, (p^(r, t)) has a quadratic 
infrared divergence. 

In our formulation of the active gel problem the diffusion coefficients D^x and Dyy for 
motions along the x- and y-directions can be calculated from the auto correlation functions 
(vj(q, oo)vj{—q^, —^)) of the velocity field v: In general Dij = | / ■^^.(Q; ^ = 0)^j(~c[, -co = 
0)). The equations of motion of Vx and Vy in the incompressible limit are given by 

Vx = -i{l'QPxxqyK.q]_ + l'2PxyqxK.q]_ " ^O^fiPxxqy " C2^fiPxyqx)Py + Pxjfj, 

.2 , D „ ,.„2 t A,.D „ t A,.D „ ^„ , D ^ (94) 



Vy = -i{l^2Pyyqxf^q± + l^0PxyqyHq± " ^2^HPyyqX " ^0AnP^yqy)Py + Pyjf. 
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where we have added thermal noise fi which is a zero-mean Gaussian noise with variances 



{fi{q±,(^)fj{-q±, -w)) = 2DkBT5i 



I]} 



(95) 



and Pij = {6ij - 



n2 . 



is the transverse projection operator. Correlators {vi{c[,uj)vj{—q, —uj)) 
can be calculated from Eqs. (1941) in a straightforward way. Evidently, off-diagonal elements 
Dxy = Dyx are zero. Coefficients D^x and Dyy have parts which depend on activity A/i and 
have infra-red divergence. Ignoring the finite parts, we obtain 



1 
2 

1 



(27r)2 



(2vr) 



X 



ql 



{vx{q±,oj = 0)vx{-q±,uj = 0)) 

[jyoPxxqyKql + V2PxyqxKq]_ + ^oAflPxxqy + ^2^fiPxyqx] ' 

A(a; = 0)r^ 



— V ^[vq sin^ 6 — U2 cos^ 6)^ 
7 D 



(96) 



where 



A(g±,a;) 



2 



^2qx) 



'M2qx + ^o^oq 



A/i 



>o6 + '^2^0) qlql + iUqx- 



(97) 



We find that if 7^ 0, then D^x has no infrared divergence. For = there are con- 
tributions to Dxx which are logarithmically infrared divergent. Further as Afi — )■ A/ic, the 
critical Afi for threshold of linear instability, Dxx diverges as (A/ic — A/i)"^. In contrast 
in equilibrium, for a passive system, Dxx is finite. The diverging contribution in the active 
system stems from the explicit dependence of Vx on py and the fact that Pj,-correlations 
are long-ranged in the plane. Moreover, using relations between the different coefficients 
as shown in Sec. IIII Al we can draw further conclusions about the diffusion coefficient Dxx 
(or, Dyy) of a tracer particle. The divergence of Dxx, now relates the critical thickness Wc 
with the critical value of the activity Afic as Wc ~ lJ_\/Afi. Thus as the A/i increases, the 
corresponding critical thickness decreases (see Ref. |32| for more details). This allows one 
to experimentally test our results in which a thickness less than Lc may be used and D may 
be measured as W ^ Wc-. Writing W = Wc - 6 with 6/Wc < 1, we obtain Dxx ~ 1/5^ 
in the limit 5 — )■ 0, when the active contribution dominates. Diffusion coefficient Dyy shows 
similar dependence on A/i or 6, although the numerical coefficient is different. 
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B. Dynamics for System II: Incompressible limit 



We begin with the constitutive relations given by Eqs. fH6|) and Wl} . As in System I (see 
Sec. IV A] above), we use the Franck free energy to calculate h. We work in the limit of equal 
Franck's constants. We have h±i = nV^Pi. We are interested in the incompressible limit 
of the dynamics and we eliminate pressure 11 using ■ v = 0. Pressure may be expressed 
as n = -^[{^2 + z^3)Vi(V± ■ hx) + (6 + 6)AAiVi(V± ■ p)]. Substituting for U, Eq. gT]) 



becomes 



Pi 



did J 



/i,- + ^ A;uP,,V"p,- + (98) 



where Pij = 5ij ^ is the transverse projection operator which appears due to the 

divergence-free condition of the velocity fields. We have added a thermal noise term gi 
which is a zero-mean Gaussian distributed noise with a variance given by 



{g,{r,t)g,{0,0))=2KBT 



^ + ^ VI P., 



5(r)5(t), (99) 



7 D 

where 6{r) is the 2d Dirac (5-function. This choice for the variances of ipi ensures that the 
FDT is held true for A/i = 0. We continue to use the same thermal noises for non-zero but 
small A/i. 

In order to study linear instabilities, we ignore the noise and separate the individual 
components of pi in fourier space and defining Pi(q, t) ~ pi{t = 0) exp(At), we find 

A = -M. _ M _ + -jkpi,i (100) 

The eigenvalues A have no angle dependence. Therefore, the 2d rotational symmetry is 
maintained. Both the eigenvalues are stable at the longest wavelengths (smallest q). Further 
one of the eigenvalues is independent of A/i and negative at all values of g, and hence, stable 
for all wavenumbers. The other one depends on A/t and may change sign (thus leading to 
instability) when z/3^2A/i < for intermediate range wavenumbers given by qWv■i^2^^^\/ D > 
^ and the crossover is determined by q'j_\h'3C,2Afi\/D = ^ yielding a crossover wavevector g^i- 
For high enough wavenumbers (dominated by the q^ term) the system is again stable and 
asymptotically matches with the equilibrium results. The crossover to stability determined 
by the q^ term is determined by the condition q\nv1/D = \i'3^2'^fJ'/ D\ defining another 
crossover scale qc2- Unlike System I, the eigenvalues have no angle dependence, and as a 
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result, in the unstable case the growth rate is same in all directions in the 2d plane. This is 
a consequence of the 2d rotational symmetry. A schematic plot of the unstable eigenvalue 
A versus wavevector q is shown in Fig. O Further, it is evident from the condition of the 
instability that for sufficiently small Afi the intermediate band of unstable wavenumbers 
may vanish. The specific condition of threshold of instability is given by qd = qc2- 



1 






q 



FIG. 9. (Colour online) A schematic plot of the eigenvalue A of the unstable mode versus wavevector 
q for System II. 

In the linearly stable case, the correlation functions of px and Py in the incompressible 
limit may be calculated exactly in a straight forward way starting from stochastically driven 
Eqs. of motion of p^ and Py (!98|l . However, we do not present the detailed but algebraically 
rather intensive calculation here. Instead we provide arguments about the correlations in 
the scaling level. Fields p^ and Py are coupled and they may be expressed in terms of 
eigenmodes whose time evolutions are independent of each other and are governed by the 
two eigenvalues of the matrix M calculated above. As we find there, one of the eigenvalues 
is always stable and independent of A/i and the other one is A/i-dependent and may be 
unstable depending upon the the value and sign of AyU. Although one of the eigenvalues is 
stable, since p^ and Py are linear functions of the eigenmodes, both of them will be affected 
by the instability, if any, of the A/i-dependent eigenvalue. Linear relations between px, Py 
and the eigenmodes ensure that the nature of divergence is same as that of the correlation 
of the unstable eigenmode. In particular there are no instabilities at the longest wavelength. 
Since the eigenvalue which depends upon Afi may change sign (thus leading to instability) 
for uo^oAfi < for intermediate range wavenumbers given by g2N^ > ^, the correlation 
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function for the corresponding eigenmode and hence the correlation functions of both px and 
Py, will show divergence when the above condition is satisfied. In the compressible limit, 
the velocity fields and Vy depend linearly on and Py] hence velocity auto-correlators 
also will show divergences for wavevectors in the intermediate range satisfying the above 
mentioned conditions. This then means that the diffusion coefficient of a tagged particle 
in the sample, being proportional to spatial integrals over appropriate velocity correlation 
functions, will have divergences once the above instability conditions are met. 

C. Dynamics for System III 

So far in the above we have set up the Eqs. of motion and analysed their instabilities 
when there are no external forces and the Eqs. for the velocity fields have a conservation 
law form. We now consider the dynamics for System III, which is the case when there are 
external forces, and, as a result, v no longer has a conservation law form. This will have 
important consequences on the ensuing dynamics as we will see below. We begin with the 
constitutive relations flM|) . flS^ and flFTl) . We consider the compressible hmit in which case 
the conservation equation is given by 



p = -V_L ■ V = D^dlU + Dydyll - uodxdyh±y - U2dyh±y - ^Q/\ixda,dyPy - ^2^fxdypll01) 

Using the thermodynamic relations 11 = (xp — aVj, ■ p) and h±y = {nV^Py + adyp) as in 
System I we get 

p = -{D^ql + Dyql)xp + ivoaqlq^p + z/gag^p + {iaD^qlqy + iaDyq^ - nu^qlqy 



where we have added zero-mean Gaussian distributed thermal noises Q with a variance 
(Ci('?,w)Cj(-?, -w)) = "^Ddij and 



i^^oqxQy + iK.i^2qyq^ + ^oA/iq^qy - i^2Afiqy)py - iq^f^ - iqyfy. 



(102) 




aqy 



7 



(103) 



where we have added a zero-mean Gaussian thermal noise fp with a variance 




(104) 
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Equations f ll02p and fllOSp have the same symmetry (polar) as those in Ref. [20]. However, 
Eqs. fll02p and fllOSp differ with the corresponding equations in Ref. j20| in details. Assum- 
ing Py(t), pit) ~ exp(At), A may be calculated easily in a straight forward way. From the 
expressions of the eigenfrequencies (we do not show the explicit forms here) we find that 
there are dissipative (or instabilities) as well as underdamped propagating modes just like 
in System / depending upon the sign of the coefficients. In the large wavevector limit, both 
the eigenvalues are stable. System III is, therefore, characterised by the presence of generic 
polar order in the ordered state, large density and polarisation fluctuations together with 



moving instabilities. In this context we refer to Ref. [26|, where the authors experimen- 
tally and numerically studied the dynamics of actin filaments in a motility assay. For high 
enough density they found a polar ordered state with larger density fluctuations coupled 
with propagating waves. We believe our formulation here is a promising starting point for 
physical understanding of actin dynamics in a motility assay. Effects of nonlinear terms are 
likely to be important in understanding the density-dependent phase transition discussed in 



Ref. [261 



When the system is linearly stable, as in Sec. IV Al we now examine the properties of rele- 
vant correlation functions in the ensuing nonequilibrium steady state. Again as in Sec. IV A[ 
such situations arise when conditions for linear instability are not met. For a compressible 
system, the density auto correlator shows novel behaviour, similar to but richer than for 
the case of System I. We obtain the equal-time density auto correlation function S{q,t) = 
/r^gS(q,u;) ^ {piq^,uj)p{-q^,-u)) as ^ (p(q, t)p(-q, t)) = ^"^^(q,^;) ~ l/q\ 
The noticeable feature is that density fluctuations are enormous as g — )■ 0, diverging as 
This would be equivalent in an equilibrium system to having a compressibility which 
diverges as l/q"^. In real space terms, this would mean compressibility x diverges as 
where L is the linear dimension of the system. As we discussed before, irrespective of the 
details of a system or any connection to a response function, the rms number fluctuations 



^y {SN'^) in an area A (we are considering a 2d system) scales as S{q — )■ 0)A. Thus for our 
system System III we have 



y/{SPP) oc y/UI oc oc iV, (105) 

where, we have used N oc A oc L"^. Thus we find giant number fluctuations in the system 
20|. 
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1. Incompressible limit 



In this limit we take V ■ v = 0. Here the only slow variable is Py whose Eq. of motion we 

set up below. The pressure can be eliminated by using the incompressibility condition and 

we get n = -^-^^^[^Ql^nq^qyPy - i^2^nqyPy - KUoq^qyq'^Py + iu2Kqyq^Py]. Using this value 

of n we write Eq. (157|) as 

nq^Py KU2q\D xcos'^ 6 — KP^q^DySm^ 6 VQ^Qq^DySm.^ 6 

^ 7 Dy{D.^ cos^ e + Dy sin^e) ^ ^ D^{D.^ cos^ 6 + Dy sin^ of^ 

- ' D^cos^ e + By sin^ 9 + '^''^ ^'''^ 

Evidently, the solution for Eq. fll06p admits underdamped waves with an anisotropic 
wavespeed C{q,9): 

gcos^sin2 0(z/2^o - i^o6)A/i 

Two situation may arise: (i) when the contribution to the propagating mode from the self- 
advecting term has the same sign as the one coming from the active terms (top figure in 
Fig. [TOl) or (ii) when they have different signs (bottom figure). In general, the wavespeed is 
anisotropic. Thus the wavespeed is zero at 9 = 7t/2 and nonzero elsewhere in this range of 
angles with a value that depends on 9 explicitly. Further, by switching the signature of the 
product (z/2^0 ~ ^0^2)'^^!' the direction of propagation can be reversed. Writing the solution 
as Py{t) ~ exp[A(g, 9)t] we find for the eigenvalue A as a function of q and 9 

^ ~ 7 Dy{D^cos'^ 9 + Dy sin^ 9) ^ D^{D^ cos^ 9 + Dy sin^ 9) 

.qcos9sm'^9{u2^o- Uo^2).. , . .,„„v 

- ' cos^ 9 + By sin^ 9 + '^''^ ^'''^ 

In general, the real part of A depends on q and 9, reflecting anisotropic growth or decay. 

Depending upon the sign of the A/x-dependent terms, there are instabilities at O(g^) when 



uo^oDyAfi 



Dx{Dx cos^ 9 + Dy sin 



2 



> 0, 



and — sm > — z , 109 

^ D^{D^ cos^ 9 + Dy sin^ 9)^ '7 Dy{D^ cos^ 9 + Dy sin^ 9)^ ^ ' 

where the equality sign determines the threshold of instability characterised by a critical 
A/ic, such that when A/i exceeds A/Xc and the above inequality is met, the system becomes 
unstable. Since the phenomenological coefficients i^o and 1^2 are in principle independent, 
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FIG. 10. A schematic polar plot of the wavespeed C{q,9) as a function of 9 for a given q. The 
anisotropic nature of the wavespeed is visible. 



such instabiUties and underdamped waves can occur independently or together. The finite 
wavevector instabihty of 0{q^) disappears at 6' = 0. At a general angle, instability may be 
present, depending upon the signs of the relevant coefficients. Therefore, as in System I, 
complicated patterns are likely to emerge out of these instabilities. The schematic form of 
the g-dependences of the unstable eigenvalue for an arbitrary 6 is similar to Fig. |6l There 
is a ^-dependent qmax at which A is maximum; the maximum instability occurs as this 
wavevector. We find 



± 



^oAfi i^2^x COS 6 Dx{Dx cos 9 + Dy sin 



(110) 



2h'oK "^i^lDy sin^ 9 2^vlDy sin'' 9 

Schematic plots of q^ax versus 9 for two different values of A/i are shown in Fig. [TTl We 
clearly see (i) qmax depends strongly on 6*, an illustration of the ensuing anisotropic pattern 
and (ii) for larger Ayu(A/i// > Afij), qmax is larger, implying that as A/x rises patterns 
become denser in the real space. The differences in the dependences of qmax on 9 indicate 
the differences in the generated patterns between System I and System III. Thus, in the real 
space a direction-dependent anisotropic pattern will emerge with a periodicity ~ ^ I qmax- A 
contour plot of eigenvalue A as a function of wavector q and polar angle 9 is shown in Fig. 
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FIG. 11. (Colour online) A schematic plot of Qmax versus 9 for some values of the parameters in 
System ///for two values of A/i : A/u// > Afij (see text). 

[T2] which depicts the stable (negative) and unstable (positive) regions. 




FIG. 12. (Colour online) A schematic 3d contour plot of the eigenvalue A as a function of wavevector 
q and polar angle 9 for some specific parameter values. Both positive and negative values of A are 
visible (see text). 

The instabihty that occurs here is for positive values of vq^q^^. Assuming positive vq, 
this would mean a constant positive active contribution ^o^Ai (using p^^ = 1) to Vx- In 
other words, when the active contribution to the velocity is parallel to the direction of the 
macroscopic order (in this case the x-direction) there is instability; when it is anti-parallel 
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to the latter, it stabilises. In terms of the relations between the different coefficients given 
in Sec. IlllCl we can get the generalised instability condition as 

where Wc is the critical thickness in the z direction at which instability sets in. It can have 
different values in different directions depending upon the values of 6'. At 6 = ^, the critical 
thickness is given by 

at which instability sets in. 

To calculate the correlation function in the incompressible limit we should add a gaussian- 
distributed thermal noise to Eq. fll06p such that it maintains FDT at AyU = 



■ 1 i4D^ ql vlDy q^y 



(113) 



.7 Dy D^ql + Dyql D^ql + Dyql 

Assuming small departures from equilibrium we continue to use the same thermal noise even 
when A/i ^ 0. This leads us to auto correlation function (p^(a;,t)) 

q±dq± dO 



{plix,t)) = kBT J 



27r 27r 



1 i^^Dx cos^ 9 ^iDy sin^ 9q\ 

7 Dy {D^ cos^ + Dy sin^ 6) D^ {D^ cos^ 9 + Dy sin^ i 



, 2/1 ^'l-Ox cos^ 6 ^iDy sin^ 6q\ 

^'^^^7 ^ Dy{D^ cos'^ 9 + Dy sin^e) ^ D^{D^ cos^ 9 + DySin^ 9)^ 
uo^oDySin^9ql 

^ D^{D^ cos^ 9 + Dysm^ 9)^ ^ ' 



Now if 



1 T^2^x cos^ 9 UQ^oDy sin^ 9q'j 



^"^^^7 ^ Dy{D, cos2 9 + Dy sin2 9) ^ D,{D, cos^ 9 + Dy sin^ 9)^^^' 
for all 9, the system is stable and {py{x,t)) has a logarithmic infrared divergence. 

The diffusion coefficient can be calculated from the velocity autocorrelation function. 
The equation for v in the incompressible limit becomes 

Vx = [iiiol^HqyPxx - t^Voqyq^Pxx) - K.l'2Pxyq]_ + C2AfiPxy]Py + Pxjfj, 

Vy = [iiCoAnqyP^y " KUQqyqlP^y) + ^2A(J,Pyy - l^2fiqlPyy]Py + Pyjfj, (116) 



where we have added a thermal noise fi with correlations given by Eq. ( ]95l) .and Pi 



{Sij — ^T^) is the transverse projection operator. We find that the diffusion coefficient D^ 
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has infra-red divergent parts which depend on A/x. Considering only the most infra-red 
divergent part, we obtain 



1 
2 



(27r)2 



{v^{q±,u = 0)v^{-q±,u = 0)) 



x[{ 



A/ig_L cos 6 sin^ 6'(z^2^o - ^^06 



1 u^^os^O ulq^sui^O 
7 ^ D ^ D 



A/iz/Q^o sin"^ ggj 



7 



z/q sin^ ^gj 



D 



D 



(117) 



Therefore, the diffusion coefficient diverges as InL, where L is the lateral system size. Thus 
we find diverging diffusion coefficient in generic polar systems, in contrast to System I, where 
such divergences are found only for nematic order. No such divergences exist in System I 
with polar order. This is directly related to the presence of the active external forces. 
Nonlinear terms, not considered here, are likely to modify the L-dependences [20]. Thus, 
Dxx diverges as where L is the linear size of the system (in the XF-plane). Further, 
near the threshold of linear instability, D^x ~ (/^c — A*)"^ where /ic the critical activity at 
the onset of linear instability in System III. 



VI. SUMMARY AND OUTLOOK 



The aim of this article is to set up a stochastically driven generalised hydrodynamic 
theory of thin confined active gels as a 2d generic model for cortical actin in eukaryotic 
cells. The presence of a substrate (the bulk of the cell) breaks Galilean invariance of the 
system. We show how effective 2d descriptions emerge for thin confined 3(i active gel systems 
by integrating the thickness dependences. This allow us to write down linear constitutive 
relations between them, subject to spatial symmetries and conservation laws. These include 
terms linear in A/x, representing small deviations from equilibrium. We show that the pairs 
{vi,Pi) and {Fi = —dill,h±i) are good candidates for thermodynamic fluxes and forces in 
the present problem. In particular, we consider three different cases - (i) in-plane nematic 
or polar order without any external force, (ii) polar order normal to the plane without any 
external force and (iii) in-plane polar order with external forces. The ensuing equations 
of motion of the relevant slow variables display linear instabilities at finite wavevectors 
for certain values and signs of the activity. In the next step, following the principles of 
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Onsager reciprocity theorem, we add thermal noises such that the FDT is maintained in 
the absence of any activity. For finite activity, the FDT is manifestly broken as is expected. 
The noisy equations of motion may display linear instability, depending upon the sign and 
magnitude of Afi. For the stable case, we use them to calculate the correlation functions of 
the appropriate variables in the nonequilibrium steady state. In addition, we calculate the 
diffusion coefficients of tagged particles for System I and III, and show that they diverge 
for a given strength of the activity. In the case when there are no linear instabilities, we 
show that the equal-time density autocorrelation function show giant fiuctuations under 



generic conditions. Our descriptions are expected to be relevant for rigid biofilms 
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of bacteria colonies on surfaces, e.g., Pseudomonus aeruginusa 3J| and microtubule arrays 



in eukaryotic cells and actin motility assays. Equations similar_to t 

bo 



lose presented here have 
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been discussed elsewhere in various contexts, see, e.g., Ref. 
here complement the existing work. Controlled cell biology experiments are needed to test 
the predictions from our models. Numerical simulations may be helpful in this regard. Giant 
:luctuations in active particle systems have already been observed in numerical experiments 
36[. We look forward to further detailed numerical work using these equations with more 
realistic situations in the future. 

Our 2d equations of motion are linear as we have neglected the nonlinear terms. In a more 
complete theoretical framework, nonlinear terms should be included. They are expected to 
be responsible in arresting the linear instabilities and ensuring transition to a non-trivial 
steady state (may be a nonuniform or a fiowing steady state). The statistical properties 



in such steady states are expected to be much richer and complex in nature 



37l |. Even in 



the case where there is no instability, nonlinear terms may modify the scaling properties 



of the correlation functions obtained at the linear level [20|, ISSj. Such issues are intimately 
connected to the general question of the nature of order in a 2d driven nonequilibrium system 
with broken continuous symmetries. Further work is in progress in this direction 27 1. 
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VIII. APPENDIX A 

Here we derive the effective 2d constitutive equations for the nematic and polar sample 
in System I directly by using symmetry arguments in 2d. As discussed above, the 2d sample 
is invariant under y — > —y as nothing distinguishes y from —y direction. However, due to 
macroscopic ordering along the x-direction (p^ = 1) there is no symmetry under x ^ —x. 
We first consider nematic order and hence p ^ — p symmetry is present. In the absence 
of any external forces, wc have v oc V(-), such that v(q = 0,t) = where q is a Fourier 
wavevector. Using these arguments we can write down the constitutive relations for Vx and 
Vy as linear functions of the relevant fluxes F^, Fy, and hy. In order to satisfy the required 
symmetry dictated tensorial nature of the constitutive relations, the Onsager matrix is to 
be constructed out of p (with — 1 here by choice) or dx, dy. In general, the constitutive 
relations for Vi have the form 

Vi = DFi + uodj{pihj) + iy2dj{pjhi) + Vdi{p ■ h). (118) 

where vq, ^2, V2 are dissipative coefficients coupling the flow with the conjugate orienting field 
h^. Next we include a nonequilibrium drive, linear in A/i. Microscopically, the conversion of 
ATP to ADP acts a chemical fuel creating an active stress (or active force) which acts on the 
fiow field over and above the equilibrium generalised forces. In addition to the equilibrium 

terms, for an active gel, there will be active contributions to f j. The most relevant form of 
such a contribution to Vi is ~ Afidj{piPj), where A^piPj is an active stress. Therefore the 
equations for Vx and Vy take the form (up to the lowest order of gradients and linear in A/i) 

Vx = DFx + i^ody{pxh±y) + iyxdx{Pxh±x) + i^2dy{pxh±y) + ^oAij,dy{pxPy), (119) 

Vy = DFy + U2dx{Pxh±y) + Vdy{pxh±x) + ^odx{Pyh±x) + ^2^lldx{PxPy): (120) 

where we have retained up the the linear order in py, Fj = — Vjll, Ux — i^o + J^2 + 1^2, 
^0 and ^2 are the coupling constants couphng activity Afi to the fiow fields. We have 
retained terms up to 0{py) in smallness in the active terms, and hence, possible active 
contributions like A^dx{p^) to Vx and A^dy{p'^) to Vy are ignored to this order. Further, 
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an active contribution of the form Afidxip^) yields zero to leading order in smallness since 
Pa; = 1 to that order. Unsurprisingly, active terms in Eqs. flll9tll20|) respect the spatial 
symmetries of the constitutive relation flllSp in equilibrium. We have used same coefficients 
^0) ^2 as in Sec. IIII Al in anticipation of obtaining equations identical to Eqs. fl34ti36p 
here. The above equations (11191) and fll20p generalise the Darcy's law for a simple fluid to 
an active (nonequilibrium) gel. 

In general in equilibrium the dynamics of p^-, Py have a term representing relaxational 
dynamics towards local equilibrium and terms which couple it to v. In this case the equation 
for Pi is 

^ = ^ + Ux,y), (121) 

where 7 is the rotational viscosity. The explicit form of the vector-valued function ipi can be 
found out using the following arguments in equilibrium (AyU = 0): Since in Eqs. f lllQp and 
f ll20p for V we have terms proportional to Vh±x and Vh±y, we expect on symmetry ground 
(Onsager principle) only gradients of v will appear in the p^ and Py equations. The signs 
and the coefficients are to be chosen in such a way that the Onsagers reciprocity theorem 
holds in equilibrium. Thus Eq. f ll2ip generally takes the form (A/i = 0), with i = x,y 

^ = ^ - aipa^d^v^ - a3PxdyVy, , (122) 
- oiop^dyVx - a2PxdxVy, (123) 



dpy _ h±y 



dt 7 

where the coefficients ai, a^, ao and 02 are to be determined using the Onsager symmetry 
(reciprocity theorem; see below). When A/i 7^ 0, explicit dependence of Vi on Afi generates 
activity- dependent terms in the equations for pi. One may further add an active term of the 
form XiAfipx in Eq. 11221 to obtain an explicit equation for p^ by substituting for Vx and Vy 
from Eqs. ffTT9|l and ffT20|) . 

dpx hj 



Lx 



Oil [yxPxdl{Pxh±x) + i^QPxdxdy{pxh^y)\ - asPx [^2dxdy{pxh±y) + iydy{pxh_Lx 



dt 7 

- aiDp^d^F^ - asDp^dyFy + XiAfip^ - ai^op^j^A/jd^dyip^Py) - ag^aPx A/i9^ (pa;P?()L24) 

However, we do not add any active term proportional to A^py directly in the Py equation, 
since Py is the transverse fluctuation and is a broken symmetry variable 3^. Substituting 
Eqs. (Iligp and (I120p in Eq. (11230 we get the equation for orientational field to the linear 
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order in gradients as 



dt 



h 



7 



±y) 



- aoDp^dyF^ - a2Dp^d^Fy 



-ao^o^fip,,dy{pxPy) - a2^2^fiPxdl{pxPy) - Px[oioJ^x + a2J^]dxdy{pxh±cc)- (125) 



We use Onsager's reciprocity relation in equilibrium (A/i = 0) to obtain 



1^0 1^2 
«0 = "2 = J^, Oil 



+ + V 

D ' 



^3 = — . 

^ D 



(126) 



In a renormalised theory coefficients and a2 should generally acquire a A/i- dependence. 
However, restricting ourselves within the scope of a linear flux-force relationship we continue 
to use the same ao, cxi-, ol2 and as as in Eqs. f ll26p for small A/i which yield the correct 
equilibrium limit. The active terms affect h^^ which may be determined as before by setting 

= 1 in Eq. fll24p to the leading order in py. We obtain 



-7 



AiA/i + -^d^dyh^y + 



^X^y^J-lf ^T.Qcr.Frp U^,C)q,Fr,i 



^y^y-^ y 



(127) 



D ^ ' D 

Substituting h±x in Eqs. (11 19p and fll20p . we note that the only effect of h_!^x is to shift the 
coefficients and .^2 of the active terms by amount Ai. Thus the nonequilibrium version of 
the Darcy's law in this case become 



Vx 



= DF^ + V0dy{p^h^y) + ^oAfldyip^Py), 

= DFy + y2dx{Pxh±y) + ^2AfJ^d^{pxPy), 



(128) 
(129) 



where coefficients ^0 and ^2 are to be understood as effective coefficients after absorbing the 
contributions coming from elimination of h±x- One may, however, add new relevant active 
terms bilinear in p and linear in V and A/i in Eq. ( ]123p : such terms however break the 
nematic symmetry considered here. We consider the effects of such terms separately below. 
Finally, putting these values of ao and a2 and eliminating h±x, Eq- (11250 becomes 



dpy 
dt 



h 



±y_ 

D 



7 



\pxdl{pxh 



±y) 



PxdliPxh 



±y) 



- i^oPxdyF^ - V2PxdxFy 



1^2 



—^lAfip^d Jp^Py) - —p^^2Afid^{p^Py). 



It is evident that the Eqs. (I128p . fll29p and fll30p are invariant under 

{i)y^~y (no distinction between?/ and — ?/), 
(m)p — )■ — p (nematic symmetry). 



(130) 



;i3i) 
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Further note that the coefficients D, z/q, 1^2 and 7 form the Onsager coefficient matrix which 
is symmetric, as is expected. The signatures of these coefficients are determined by the 
positivity of dissipation. Substituting Eqs. f ll281ll29|) and f ll30p in Eq. ([T]) and demanding 
that each term makes a positive contribution to dissipation separately, we find D > 0, 7 > 
0, Uq > and 1^2 > 0. 

The above equations (I128p . (I129P and fllSOp exhibit the nematic order. In order to 
introduce polar order this inversion symmetry must be broken. Hence we add a term ^a{p ■ 
V_l)p to the equation for py. Since = 1, the only term, which survives upon linearising 
about the chosen reference state, is ^a^- Here, is an active coefficient, proportional to 
A/i. Hence the equation for py having polar order is 

- ^oPxdyFx - V2PxdxFy 
-^^O^f^PxdyiPxPy) - ^^2^f^PxdliPxPy) + ^APx^- (132) 

Note that Eq. (11320 breaks the invariance under p — ?► — p and (separately) under r — )■ — r. 
The lack of the former makes the system polar; the lack of the latter leads to propagating 
wave- like excitations, as we will see below. Finally, using p^. = 1 we obtain Eqs. f l34H36|) . 



Py 



- T^P^d'ipxh^y) - 



7 



D' 



D 



PxdliPxh 



±y) 



IX. APPENDIX B 

Here we derive the 2d equations directly for a sample with external forces. In this sample 
the total generalised forces, due to the presence of external fores, is not zero and must be 
proportional to either h±[q = 0,t) or p(q = 0, t)A/i (q is a fourier wavevector). As the 
system is polar we do not have p — ?■ — p symmetry. We further consider planar alignment, 
i.e., in this case we have (p) = po^x = Bx 7^ in the plane. Hence, there is no x — )■ —x 
symmetry. However, the system is still invariant under y — )■ —y. We continue to impose 

= 1. In addition we consider departure from equilibrium to the linear order in Afi. We 
consider a situation where the system experiences external forces (e.g., forces imparted to 
the actin ffiaments by the immobile molecular motors grafted on the confining substrate in 
actin motility assays) which, on symmetry ground, should be proportional to hi or A/xpj. 
From the invariance properties discussed above, we can write the constitutive relations for 
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the velocity fields such that the Onsager reciprocity principle holds in equilibrium (A/x = 0) 

DF^ + i^ody{p^h±y) + U2h±^ + ^oAfxdyip^Py) + ^oAfidyip^Py), (133) 
.y - DFy + U2h±y + ^2^fiPy, . (134) 



where we have retained terms with leading order spatial gradients, up to linear order in 
Py and active terms up to linear order in A/i having same symmetries as the equilibrium 
terms. Coefficient D > is the inverse friction coefficient, and z/q, 1^2 and Uy are dissipative 
cross-coupling coefficients, coupling flows with the local orientation, ,^0 and .^2 are coefficients 
of the active terms. 

Equilibrium dynamics of the orientational field comes from its relaxation to local equi- 
librium and its coupling to flows. Hence the general form of the orientational field can be 
written as (to leading order in gradients) 



Px 



Py 



h 



7 
7 



+ a^Vx + ayP^dyVy, 
+ a^Vy + a2PxdyVx. 



(135) 
(136) 



where ax, ao, a2 and ay are coupling constants coupling the orientation field px, Py to the 
flows, all of which can be determined by using the Onsager reciprocity theorem in equilib- 
rium, i.e., A/i = 0. When the active terms in Eqs. fll33p and fll34p are non-zero, explicit 
dependences of Vx and Vy generate A/i-dependences of px and Py. In addition we add an 
active term XiAfipx in the p^-equation. However, a similar term does not exist in the dy- 
namical equation of Py, as to the lowest order py is the transverse mode. Put together 
everything, we get (to the leading order in gradients) 
. hi 



Px 



7 



+ ax[DFx + iy2h±x + Jyody{pxh±y)] + ayPxdy[DFy + iy2h±y + iyy{pxh±x)] 



+ XiAfipx + ax^oAndyipxPy) + ay^2^fipxdyPy, 
Py = + "o'^2 - a2i^ody]hj_y + a^DFy + a2DdyFx + a2i^2Pxdyh±x 
+ a2^oA^pxdyPy. 



(137) 



(138) 



The condition px = 1 yields h±x, which upon substituting in Eqs. fll33p and fll34p give (for 
A/i = 0): 



D 



1 - 



h + <^xV2 



Fx + z^o 



1 - 



^ + 01x^2 



dy{pxh 



Vy = DFy + U2h±y. 



(139) 
(140) 
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Thus the effective friction coefficients are no longer isotropic. Further, Onsager reciprocity 
theorem requires 

1^2 ^0 ^2 1^0 

ao = j^, a2 = - — , = —, ay = - — . (141) 
Defining = D'^/{D + z/|7), Dy = D and u'q = VqD/{1 + 7Z/|), we write 

v^, = D^F^ + UQdy{p^h±y), (142) 

Vy = DyFy + 1^2^^ ±y , (l43) 

Evidently, the fiction coefficients are now anisotropic and < Dy. This is in contrast to 
the nematic sample or System I, which we considered above, where fiction coefficients are 
isotropic. For the purpose of brevity, we represent z/q by uq. As in Appendix A, we continue 
to use the same equilibrium values for ao, 0:2 and ax even when Afi ^ 0. Then using their 
values together with h±x and setting p^; = 1 we then obtain the constitutive relations ([541 - 
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